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Abstract. New results from the K2K (KEK to Kamioka) long-baseline neutrino- 
oscillation experiment is presented 



1. Introduction 

> 

The K2K experiment [T], |2], |3], f|, g is the first long-baseline neutrino-oscillation 
experiment with a distance of hundreds of km and using an accelerator-based neutrino 
beam. The nominal sensitive region in the neutrino-oscillation parameters is Am 2 > 
3 x 10~ 3 eV 2 . This covers the parameter region suggested by the atmospheric neutrino 
anomaly observed by several underground experiments, || [7], |8| and confirmed by Super- 
Kamiokande (SK)||. 

An overview of the K2K experiment is as follows. An almost pure wide-band z/ M 
beam generated in the KEK 12-GeV Proton Synchrotron (PS) is detected in the SK 
detector at a distance of 250 km. Two beam monitors are located along the beam 
line. One is a pion monitor (PIMON), which is a gas Cherenkov detector installed 
downstream of the second magnetic Horn. The other is a muon monitor (MUMON), 
which is an ionization chamber and a silicon pad array behind the beam dump. In 
addition, two different types of front detectors (FDs) are located at the KEK site. One 
is a lkt water Cherenkov detector (1KT), which is a miniature of the SK detector. The 
other is a so-called fine-grained detector (FGD), which is composed of a scintillating 
fiber tracker (SFT)fL0|, trigger counters (TRG), lead glass counters (LG) and a muon 
range detector (MRD)|T|. 

Since the neutrino beamline, the beam monitors and detectors are already described 
in previous articles|3|, f|, |5|], they are not described here. Instead, summaries of these 
detector components are given in tables. Table 1 gives the design features of the neutrino 
beamline. Table 2 gives the design and performance features of beam monitors. The 
design, purpose, and performance of the detectors are given in Table 3. A schematic 
view of the Front detectors is also shown in Fig.l. 

The K2K experiment was successfully started in early 1999, and physics data were 
recorded from June 1999 to July 2001. The total data-taking period during three years 
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Table 1. Design of the neutrino beam line. 



Proton kinetic energy 
Proton intensity 
Extraction mode 
Beam duration 
Target 

Decay tunnel 



12 GeV 

5.6 x 10 12 proton/pulse 
Fast extraction 

~ 1.1 /^second for every 2.2 second 
Aluminum (3cm0x 65cm) 
200 m 



Table 2. List of the beam monitors along the K2K neutrino bcamlinc. 
Pion monitor (PIMON) 

• Design Gas Cherenkov detector with a spherical mirror and R-C318 gas. 

Reflected light is collected by 20 Hamamatsu R5600-01Q PMTs. 

• Position Downstream of the second magnetic horn. 

• Purpose Cherenkov photons from charged pions are measured 

at several refraction indices controlled by the gas pressure. 
The momentum and angular distributions of the pion beam are obtained 
from the intensity and shape of the Cherenkov image in the focal plane. 
The neutrino energy spectrum for E v > lGeV can be calculated 
from the kinetics of the pion decay. 
Muon monitor (MUMON) 

• Design 2m x 2m ionization chamber filled with He gas 

and a silicon pad array. 

• Position Downstream of the beam dump. 

• Purpose The position of the beam center is monitored with an accuracy of 

less than 2 cm in spill on spill basis. 



Figure 1. Front detectors in the K2K experiment. 
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Table 3. Design, purpose, performance and event rate of the detectors used the K2K 

experiment. 
Front detectors 300m downstream 
(A)lkt water Cherenkov Detector 

• Design A 1/50 (total volume) miniature of the Super-Kamiokande detector; 496 tons 

of water are viewed by 680 20-inch PMTs. The fiducial mass is 25.1 ton. 

• Purpose Direct comparison with neutrino events in Super-Kamiokande 



• Performance e//i identification capability: >99%; AE e /E e : 3%/yE e (GeV) 

• Event rate ~0.005 events/spill in fiducial volume; ~1 xlO 5 events for 10 20 p.o.t. 
(B)Fine Grained Detector 

• Design Consists of a scintillating fiber tracker, trigger counters, 

lead-glass counters and a muon range detector. 

• Purpose Precise measurement of the neutrino profile and energy distribution 
(B-l) Scintillating fiber tracker (SFT) 

• Design 20-layer "sandwich" of scintillating fiber (O.7mm0) sheets 

and water in aluminum containers. 

Sensitive area: 2.4m x 2.4m; Fiducial mass: 5.94 tons 

• Function Track reconstruction of charged particles and identification of the kinematics 

of neutrino interactions. Water is used as the target material. 

• Performance Detection efficiency of each layer: >99% 

Position resolution: ~ 0.6mm 

• Event rate 0.001 events/spill, 2 xlO 4 events for 10 20 p.o.t. 
(B-2)Trigger counters (TRG) 

• Design 80 large plastic scintillators (466cm (L)xl0.7cm(W)x4.2cm(T)) 

covering upstream and downstream of the scintillating fiber tracker 

• Function Rejection of cosmic ray muons and neutrino events from lkt water 

Cherenkov detector, and measurement of the absolute event time. 

• Performance Timing resolution: ~2nsec; Position resolution: ~5cm. 

Detection threshold@center of the scintillators: ~1.7MeV. 
Detection efficiency for a penetrating charged particle: > 99%. 
(B-3)Lead glass counters (LG) 

• Design 600 lead glass counters with an acceptance of 11.3cm xl2.2 cm each 

• Function Identification of electrons from the energy deposit in the counter 



• Performance AE e /E e : ~ 10%/^ £ e (GeV). 
(B-4)Muon range detector (MRD) 

• Design 12-layer "sandwich" of ~900 drift chambers 

and iron filters (10cm-20cm thickness). 

• Function Measurement of the muon energy from the range. 

• Performance Position resolution: 2.2mm. 

Detection efficiency of each layer: ~ 99%; AE^/E^: 8 ~ 10%. 
Far detector (Super-Kamiokande) 250km downstream 

• Design 50kt huge water Cherenkov detector at about 1000m underground. 

22.5kt of the fiducial volume is viewed by 11164 20-inch <p PMTs. 

• Performance e//i identification capability: > 99 %; 



AE e /E e : 3%/ v /£ e (GeV); AEJE^. 3% 
Accuracy of the absolute event time adjustment: <0.2/zsec 
• Event rate ~0.3 events/day, ~160 events for 10 20 p.o.t. 
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was 234.8 days. The accumulated beam intensity was 47.9 x 10 18 protons on the target 
(p.o.t.), which is about 50% of the goal of the experiment, 10 20 p.o.t. 



2. Study of the Neutrino beam properties at the KEK site 

The characteristics of the neutrino beam at the KEK site were examined using FDs and 
beam monitors. In this section, the present status of analyses on (l)the neutrino beam 
direction, (2)the neutrino beam intensity and its stability, (3)the z/ e /V M ratio, (4)the 
neutrino energy spectrum, and (5)extrapolation of the neutrino flux are presented. 



2.1. Neutrino beam direction 



The neutrino beam-line was constructed using a GPS position survey p2| ; the alignment 
of the beam-line, FDs and SK is better than 0.1 mrad. The neutrino beam direction 
relative to the beam-line was measured with MUMON and MRD independently. 

The position resolution of MUMON is about 2 cm, corresponding to an angular 
resolution of 0.1 mrad. Because and muons originate in the same pion decay in the 
decay volume, the beam direction can be examined from the profile center of the 
muon beam. The time variation of the profile center is plotted in Fig. 2. The direction 
of the muon beam agrees with the beam-line within 1 mrad. 

The neutrino beam direction is also measured using neutrino interactions in 



MRD [|TT|j . From spatial constructions of the vertex positions, the center of the vertex 
positions is found to agree with the SK direction within 1 mrad. The time variation of 
the beam center, also plotted in Fig. 2, shows that the steering of the beam direction is 
stable, and is consistent with the results from the muon monitors. 

Since the angular acceptance of the SK detector from the KEK site is about 
0.2 mrad, and the energy spectrum of the neutrino beam is expected to be uniform 
within 3 mrad from the center of the beam axis, the adjustment of the neutrino beam 
direction (< 1 mrad) is sufficient. 
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Figure 2. Time variations of the beam directions measured by MUMON (left) and 
the time variations of the profile center measured by the vertex distributions in MRD 
(right). The top figures and the bottom figures are the horizontal distributions and the 
vertical distributions, respectively. The SK direction and 1 mrad off axis are shown by 
the solid and dashed lines, respectively. 
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2.2. Neutrino beam intensity and its stability 

The absolute neutrino event rates in the FDs are used to calculate the expected number 
of events in SK. The neutrino events in 1KT are reconstructed from the Cherenkov 
ring pattern by the same method as that used in the SK atmospheric neutrino analysis. 
We selected events with the reconstructed vertex being inside 25.1 tons of the fiducial 
volume, which is defined by a 2 m radius, 2 m long cylinder along the beam axis. 
Since the average number of neutrino events in the whole 1KT detector is about 
0.2 events/spill, multiple neutrino events in one spill must be taken into consideration. 
The total sum of PMT signals, which are recorded by FADC, is used for counting the 
number of neutrino interactions in a spill. By considering this correction, the number 
of neutrino interactions in the fiducial volume was obtained to be ~0.02 events per 
spill. This measurement has a 4% systematic uncertainty, which is mainly due to vertex 
resolution of the reconstructions. 

Although the absolute event numbers were also studied by SFT and MRD, the 
neutrino events in 1KT are preferable for an absolute event number analysis, because 
the possible systematic errors inherent to the detection method are the same as that of 
SK, and are canceled in the analysis. Therefore, details of the SFT and MRD analysis 
are not presented in this report. 

The stability of the beam intensity is contaneously measured from the neutrino 
event rate in MRD because of their large statistics. A detailed description of the event 



reconstruction can be found in [111. The time variation of the neutrino interactions in 



MRD is shown in Fig. 3. The neutrino beam intensity is stable within a few %. 
2.3. Ue/Ufj, ratio 

The Ve/v^ ratio of the neutrino beam at the KEK site was measured by 1KT and FGD. 
The idea of the measurements is given in ]3J. An analysis with 1KT is still under way, 
and no numerical result has been obtained yet. On the other hand, a very preliminary 
result with FGD is reported to be (1.8±0.6 + 1 q)%RI3"11, where the expectation based on 
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Figure 3. Time variation of the neutrino event numbers in 1KT. The June 99 data is 
smaller than the other periods because the current of the magnetic horn was different. 
The neutrino beam intensity is stable within the statistical errors. 
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a Monte- Carlo simulation is 1.3%. The original neutrino beam has been proved to be 
almost a pure beam. 

2.4- Neutrino energy spectrum 

The neutrino energy spectrum was examined with 1KT and FGD. The 1KT has a high 
efficiency for low-energy neutrinos below 1 GeV and a full 4tc coverage in the solid angle. 
However, the efficiency for muons above ~1 GeV is poor because most of the muons 
escape from the detector. On the contrary, although FGD has a high efficiency for 
higher momentum muons, the efficiency for low-momentum muons is relatively poor, 
and the angular acceptance is limited to be within 50 degrees from the beam axis. By 
combining the 1KT and FGD results, the entire range of the neutrino energy spectrum 
can be covered. 

In the 1KT detector, events with only one /i-like Cherenkov ring are selected. It was 
also required that the vertex position is in the 25.1 tons of the fiducial volume described 
in the previous section, and the particles are fully-contained in the inner detector. To 
determine the neutrino energy, single /z-like ring events are assumed to be the charged- 
current quasi-elastic (CCQE) interaction of muon neutrinos, v^N — > fiN'. Most of 
the neutrino energy is transfered to the muons in this interaction mode. The neutrino 
energy (E v ) can be directly calculated from the momentum (p^) and angle of the 
outgoing muons as 

m N E^ - m 2 /2 
m N - E^ + p^ cos M 

where and are the masses of the nucleon and the muon, respectively, and 
E^ = y/p^ + m 2 ^. The muon momentum and angle from the beam axis are shown 
in Fig. 4 together with the Monte Carlo expectations described below. 

In FGD, neutrino events which satisfy the following criteria are selected: (l)the 
vertex position is inside of the 5.94 ton SFT fiducial volume; (2)one track passes 
through 3 layers of scintillating fibers, TRG, LG, and is stopped in MRD. CCQE events 



1 kt: ^-momentum Distribution (Fid.25t FC 1-Ring ji-like) 
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Figure 4. Muon momentum and angle from the beam axis at the 1KT detector. The 
cross is data and the box is a Monte-Carlo simulation with the best-fit parameters. 



Present Status of the K2K Experiment 



7 



are recognized as 1-track or 2-track events because the recoiled proton track can be 
occasionally identified as another short track. If the CCQE interaction is assumed, the 
travel direction of the recoiled proton can be expected from the direction of the muon 
track by using momentum conservation. The angle between the expected direction 
and the real second track (AO) can be used to distinguish CCQE events from other 
interactions. Three categories of events are defined as: (1) 1-track events, (2)2-track 
events with AO < 25°, and (3) 2-track events with AO > 30°. The second and the 
third samples are CCQE and non-CCQE enriched event samples, respectively. The 
momentum (p M ) and angle (6^) distributions for three events samples are shown in 
Fig. 5. 

The expectations are based on GEANTfQ with a detailed description of the 
materials and magnetic fields in the target region and the decay volume. It uses 
measurements of the primary-beam intensity and profile at the target as input 
parameters. The primary proton interactions on aluminum are modeled with a 
parameterization of hadron production data[|HJ. Other hadronic interactions are treated 
by GEANT-CALOR0. 

The momentum distribution in 1KT and FGD can be used to determine the 
neutrino energy spectrum by the following procedure. The neutrino energy spectrum 
is divided into 8 energy regions, and the corresponding muon momentum distributions 
are prepared by a Monte-Carlo simulation for each detector. The weighting factor of 
each region as well as other systematic errors are assumed to be fitting parameters. By 
the comparing with the real data in 1KT and FGD, the neutrino energy distribution 



SciFi P.. distributions 



SciFi &„ distributions 



1500 




* 1 Track 








1000 




• Data 








500 




^Ctfc, □ (CCQE) 
. ' ■ . i i ' T i » l i i i 


u o 


1 2 3 ' 






(GeV/c) 



300 
200 
100 




200 






2 Track 


150 




+ th 


A0 p < 25° 


100 








50 












I j i i_ I i_ 




"0 


1 2 


3 i 








(GeV/c) 



2 Track 

A0 p > 30 { 



1 



P, (GeV/c) 




1 Track 

• Data 
□ MC 
Q (CCQE) 

i l i i_ 

75 

(deg.) 



2 Track 

A0 p < 25° 



75 

(deg.) 



200 
100 




2 Track 

A0 p > 30° 



25 



50 



75 
% (deg.) 



Figure 5. Muon momentum and angle from the beam axis at the FGD detector. The 
cross is data and the box is a Monte-Carlo simulation with the best-fit parameters. 



Present Status of the K2K Experiment 



8 



is determined. The most probable values and their possible ranges of the systematic 
errors are also obtained from the fitting, and are applied to the following SK analysis. 

The best-fit results are also shown in Fig. 5 and Fig. 6. The agreement with the 
data is excellent. 



2. 5. Extrapolation of the neutrino flux 

The production of neutrinos is calculated from the simple two-body decay of the charged 
pions. Therefore, once the kinematic distribution of charged pions after the second 
magnetic horn is known, it is possible to predict the muon neutrino spectrum at any 
distance. 

The pion momentum and angular distribution are measured by PIMON, whose 
characteristics are given in Table 2. In PIMON, intensity and shape of the Cherenkov 
image in the focal plane are measured by the PMT array. The measurement is a 
superposition of slices of the Cherenkov ring from charged pions of various momentum 
and angles. The refraction indices are controlled by the gas pressure, and seven 
independent measurements at different Cherenkov thresholds are done. From these 
measurements, the relative pion beam intensities are obtained as a function of the 
momentum and angle. To avoid contamination of the Cherenkov light from the surviving 
12-GeV primary protons, the refraction indices are adjusted to be below the Cherenkov 
threshold of the 12-GeV proton. Accordingly, an analysis for the pion momentum less 
than 2 GeV is not possible, and the measurable neutrino energy range is limited to be 
larger than 1 GeV. A more detailed description of the PIMON and the analysis method 



can be found in 117 



Fig. 6 shows the energy spectrum shape and the flux ratio in SK and FDs (so- 
called Far/Near flux ratio) obtained by the above-mentioned analysis together with the 
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Figure 6. The top figure shows the spectrum shape of the muon energy at FDs. The 
bottom figure shows the far to near flux ratio. The histograms are from the beam 
simulations. The data points are derived from the PIMON measurements. 
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results from a Monte-Carlo simulation. The agreements between this analysis and the 
simulation are excellent. Fig. 6 (bottom) can be used to calculate the expected event 
numbers in SK from the observations in the FDs. 



3. Observation in Super-Kamiokande 



To obtain beam-correlated fully contained neutrino interactions, an event selection 
similar to an atmospheric neutrino analysis was applied. For the selected events, 
the time correlation with the neutrino beam was then examined. Fig. 7 shows the 
time difference between the neutrino beam and the events obtained from atmospheric 
neutrino selection. 

Considering the neutrino beam duration (1.1/xsec) and accuracy |18|] of the absolute 
time determination (< 0.2/xsec), events within a 1.5/isec time window covering the 
neutrino beam period were selected. A total of 56 fully-contained events were found 
in 22.5kt of the fiducial volume. The number of the neutrino events and the results 



from the particle identification [19[ are summarized in Table 4. Because the expected 
atmospheric neutrino background in the fiducial volume within the neutrino beam period 
was calculated to be 1.3 x 1CT 3 events, the 56 events in the fiducial volume are a clear 
signal of neutrinos from KEK. 




A(T)(is 



Figure 7. Time correlation between the neutrino beam period and SK events which 
are selected by the standard atmospheric neutrino analysis. Events in the 1.5 //sec 
gate (—0.3 /isec ~ 1.2 /xsec) are finally selected. 



4. Oscillation analysis 

Strategies of the «-> v T oscillation study at K2K are summarized as follows. If 
Up oscillate to u T , the absolute event number of neutrino interactions in SK is reduced 
because v T interact only through neutral current interactions in the K2K neutrino energy 
range. In addition, the neutrino energy spectrum in SK should be distorted, because 
the oscillation probability is a function of the neutrino energy. 

In the following two subsections, the status of the data analyses on the absolute 
event numbers and on the distortion of the neutrino energy spectrum are presented. In 
the final subsection, an analysis of the constraints on the neutrino oscillation parameters 
is discussed. 



4-1. Absolute event numbers 

The expected event numbers in SK are calculated from the neutrino event rate in FDs, 
and an extrapolation from the FDs to SK. The numbers obtained in 1KT, SFT, and 
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Table 4. Number of neutrino events in SK. Expectations based on the event rate at 
1KT, SFT, and MRD are also shown. 



Event Category 


SK data 


(1KT) 


Expected 
(SFT) 


(MRD) 


Single ring events 


32 








(/i-like) 
(e-like) 
Multi ring events 


30(29) 
2 
24 








Total 


56 


80.1+^ 


87 5 +my 


87 4 +ri - Y 

° 1 -13.9 



MRD are used for the event rate in the FDs, as reported in 2.2. The Far/Near flux ratio 
has already been obtained from PIMON, as presented in 2.5. The expectations based on 
the data from the FDs are 80.1+JJ(1KT), 87.5™ (SFT), and 87.4+}^ (MRD). These 
results are consistent with each other. As mentioned in 2.2, we used the numbers from 
1KT as official numbers because of its small systematic errors. 

The statistical probability that the observation is equal to or smaller than 56, where 
the expectation is 80.1 is about 1%. The observation is significantly smaller than 
the expectation. 

4-2. Shape of the neutrino energy spectrum 

The neutrino energy spectrum is calculated from 29 single-ring /i-like events. (Among 
30 single-ring /i-like events listed in Table 4, one event recorded in the early stage of the 



20 




K< ec (GeV) 

Figure 8. Neutrino energy spectrum reconstructed from 29 single /i-like events 
in SK shown by points with error bars. The blue boxes are the expectations in 
the case of null oscillation. The sizes of the boxes are the systematic errors of 
the expectations. The red histogram is the expectation for the best-fit oscillation 
parameters; (Am 2 , sin 2 29) = (2.8 x l(T 3 cV 2 ,l). 
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experiment is not used in the spectrum analysis because the condition of the neutrino 
beam line is different from that of the other data-taking period. See also Fig. 3.) 

Fig. 8 shows the energy spectrum of the 29 events with expectations. These 
expectations were obtained from the neutrino energy spectrum discussed in 2.4, and 
the Far /Near flux ratio presented in 2.5. A reduction of the neutrino event in the 
0.5-1.0 GeV bin is obvious. If the absolute normalization of the flux is assumed to be 
a free parameter without any constraints, the probability that the data is explained 
by statistical fluctuation of the null oscillation is calculated to be about 15%. A more 
quantitative analysis is done by combining with the absolute event numbers. 



4-3. Constraints on the oscillation parameters 

To evaluate the absolute number of events and the spectrum shape simultaneously, an 
un-binned maximum likelihood analysis is employed. The likelihood is defined as the 
product of the normalization term and the shape term: 

£(Am 2 , sin 2 29) = C norm {Am 2 , sin 2 29) x £ shape (Am 2 , sin 2 29). (2) 

The normalization term, L norm , is a simple poisson probability of observing N b s (=56) 
for the prediction N exp (Am 2 , sin 2 29). The shape term, C s h ape , is the product of the 
probabilities for each single /i-like events: 

C shape (Am 2 } sin 2 29) = \[ Vi (E t - Am 2 , sin 2 29); (3) 

i=l 

where pi is the probability of i-th event having energy Ei at a given set of oscillation 
parameters. iVi M is the number of single /i-like events used in the analysis, and is 29. 

When the likelihood, C, is maximized, various systematic errors are taken into 
consideration. The source of the systematic errors include the spectrum shape, the 
Far/Near flux ratio, non-CCQE/CCQE ratio, the efficiency in SK, the absolute energy 
scale in SK and so on. Not only the individual systematic errors, but also the correlation 
between them are carefully taken into account. For the systematic error parameters 
which are common with FD, their probable range from the FD analysis are considered. 
More details about the systematic errors and their threatments are not discussed in this 
report, but is reported are [E0 . 



The oscillation parameter set which maximize the likelihood, C, was calculated to 

be 

(Am 2 , sin 2 29) = (2.8 x 10 _3 eV 2 , 1.0). (4) 

In this case, the expected total event number is 54.2. The expected energy spectrum 
is shown by the red line in Fig. 8. The agreements with the experimental results are 
excellent. 

The probability that the K2K results can be explained by a statistical fluctuation 
of null oscillation is calculated from the likelihood ratio of the null oscillation to the 
best-fit point. The results were obtained to be 0.4% or 0.7%, depending on the two 
independent treatments of the systematic errors in the likelihood |2(J. We conclude that 
the null oscillation probability is less than 1%. 

Finally, constraints on the oscillation parameters are studied from the likelihood 
ratio to the best fit. The allowed region in the oscillation parameter plane is shown in 
Fig. 9. The 90% C.L. allowd region on Am 2 for the full mixing, sin 2 29 = 1, is 

Am 2 = (1.5 ~ 3.9) x 10~ 3 eV 2 . (5) 

These results agree well with the most recent SK atmospheric neutrino analysis |2l], as 
shown in Fig. 9. 
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5. Summary 

The K2K long-baseline neutrino-oscillation experiment has been successfully operated 
since 1999. By the end of 2001, a total intensity of 47.9 xlO 18 protons on target was 
accumulated. A total of 56 fully-contained neutrino interactions in the 22.5kt of the 
fiducial volume of the SK detector were observed, where the expectation based on the 
data from the Front Detectors is 80. l]^' 2 . The neutrino energy spectrum obtaind from 
29 single /i-like ring events shows a clear deficit in neutrino energy between 0.5 GeV and 
1.0 GeV. The probability that the measurements at SK can be explained by a statistical 
fluctuation of the null oscillation is less than 1%. The 90% CL. allowed region on the 
oscillation parameters is Am 2 = (1.5 ~ 3.9) x 10 _3 eV 2 for sin 2 29 = 1. This result is 
consistent with the results from the SK atmospheric neutrino analysis. 
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Figure 9. The 68%, 90% and 99% C.L. allowed region on the neutrino oscillation 
parameters for <-> v T oscillation. The K2K results (solid lines) are overlaid with the 
SK atmospheric neutrino results (dashed lines). 
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